MOA 2003-BLG-37: A Bulge Jerk-Parallax Microlens Degeneracy 



B.-G. Park 1 , D.L. DePoy 2 , B.S. Gaudi 3 , A. Gould 2 , C. Han 2 ' 4 , and R.W. Pogge 2 

(The //FUN Collaboration) 
and 

F. Abe 5 , D.P Bennett 6 , LA. Bond 7 , S. Eguchi 5 , Y. Furuta 5 , J.B. Hearnshaw 8 , K. Kamiya 5 , 
P.M. Kilmartin 8 , Y. Kurata 5 , K. Masuda 5 , Y. Matsubara 5 , Y. Muraki 5 , S. Noda 9 , K. 
Okajima 5 , N.J. Rattenbury 10 , T. Sako 5 , T. Sekiguchi 5 , D.J. Sullivan 11 , T. Sumi 12 , P.J. 
Tristram 10 , T. Yanagisawa 13 , and P.C.M. Yock 10 
(The MOA Collaboration) 



ABSTRACT 



1 Korea Astronomy Observatory, 61-1, Whaam-Dong, Youseong-Gu, Daejeon 305-348, Korea; bg- 
park@boao.re.kr 

2 Department of Astronomy, The Ohio State University, 140 West 18th Avenue, Columbus, OH 43210, 
USA; depoy, gould, pogge@astronomy.ohio-state.edu 

3 Harvard-Smithsonian Center for Astrophysics, Cambridge, MA 02138, USA; sgaudi@cfa.harvard.edu 

4 Department of Physics, Institute for Basic Science Research, Chungbuk National University, Chongju 
361-763, Korea; cheongho@astroph. chungbuk. ac.kr 

5 Solar- Terrestrial Environment Laboratory, Nagoya University, Nagoya 464-8601, Japan; abe, furuta, 
kkamiya, kmasuda, kurata, muraki, okajima, sado, sako, sekiguchi, ymatsu@stelab.nagoya-u.ac.jp 

6 Department of Physics, Notre Dame University, Notre Dame, IN 46556, USA; bennett@emu.phys.nd.edu 

institute for Astronomy, University of Edinburgh, Edinburgh, EH9 3HJ, UK; iab@roe.ac.uk 

8 Department of Physics and Astronomy, University of Canterbury, Private Bag 4800, Christchurch, New 
Zealand; john. hearnshaw, pam.kilmartin@canterbury.ac.nz 

9 National Astronomical Observatory of Japan, Tokyo, Japan; sachi.t.noda@nao.ac.jp 

10 Department of Physics, University of Auckland, Auckland, New Zealand; nrat001@phy.auckland.ac.nz, 
paulonieka@hotmail.com, p.yock@auckland.ac.nz 

11 School of Chemical and Physical Sciences, Victoria University, PO Box 600, Wellington, New Zealand; 
denis . sulli van@ vuw . ac . nz 

12 Department of Astrophysical Sciences, Princeton University, Princeton NJ 08544, USA; 
sumi@astro.princeton.edu 

13 National Aerospace Laboratory, Tokyo, Japan; tyanagi@nal.go.jp 



- 2 - 



We analyze the Galactic bulge microlensing event MOA 2003-BLG-37. Al- 
though the Einstein timescale is relatively short, t E = 43 days, the lightcurve 
displays deviations consistent with parallax effects due to the Earth's accelerated 
motion. We show that the x 2 surface has four distinct local minima that are in- 
duced by the "jerk-parallax" degeneracy, with pairs of solutions having projected 
Einstein radii, f-E = 1.76 AU and 1.28 AU, respectively. This is the second event 
displaying such a degeneracy and the first toward the Galactic bulge. For both 
events, the jerk-parallax formalism accurately describes the offsets between the 
different solutions, giving hope that when extra solutions exist in future events, 
they can easily be found. However, the morphologies of the x 2 surfaces for the 
two events are quite different, implying that much remains to be understood 
about this degeneracy. 

Subject headings: gravitational lensing 



1. Introduction 

Microlens parallaxes are playing an increasingly important role in the analysis of mi- 
crolensing events. The microlens parallax, n E is a vector, whose magnitude is the ratio of 
the size of the Earth's orbit (1 AU) to the Einstein radius projected onto the plane of the 
observer (f E = AU/Ve), and whose direction is that of the lens-source relative motion. From 
simple geometrical considerations (e.g., Gould 2000), 



7r re i „ /—— 4G o mas 

** = M^ e E = ^M^ h K = — c8.1 — , (1) 

where M is the mass of the lens, 7r rc i is the lens-source relative parallax, and 9 E is the angular 
Einstein radius. Hence, if both n E and 9 E can be measured, one can determine both M and 
7r re i. While this has only been done for one microlensing to date (An et al. 2002), there 
are plans to apply this technique to several hundred events using the Space Inter ferometry 
Mission (SIM) (Gould & Salim 1999). 

Moreover, even in cases for which 9 E cannot be measured, the microlens parallax can give 
very important information about the event. When combined with the Einstein timescale, 
t E , it yields v, the lens-source relative velocity projected onto the observer plane, 

7T E £ E 
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By combining measurements of 7Te and £e with known Galactic structure parameters, Mao 
et al. (2002) and Bennett et al. (2002) were able to show that several microlensing events 
seen toward the Galactic bulge are most probably due to black holes. 

Microlens parallax measurements are possible whenever can be compared to some 
"standard ruler" in the observer plane. The two logical possibilities are observing the event 
simultaneously from two different locations and observing the event from an accelerated 
platform. Since is typically of order an AU, the most practical realizations of these ideas 
would be, respectively, to observe the event simultaneously from a satellite in solar orbit 
(Refsdal 1966) and (for sufficiently long events) to observe it from the Earth as it moves 
about the Sun (Gould 1992). 

When Refsdal (1966) first proposed obtaining microlens parallaxes by combining space- 
based and ground-based observations of microlensing events, he already noted that these 
measurements would generically be subject to an ambiguity. Using somewhat different no- 
tation, he showed that, 

AU 

tte = (7Te,||, 7r E> ±) = -t— (Ar, Au ) (3) 

"sat 

where d sat and the parallel (||) direction are defined by the magnitude and direction of the 
Earth-satellite separation vector projected onto the plane of the sky, Ar = At /t E is the 
difference in times of maximum, to, as seen from the Earth and satellite (normalized to t^), 
and Auo is the difference in impact parameters. While Ar (and so 7Te,||) is unambiguously 
defined, 

Au = ±|uo 1 e±«o, 8 at|, ( 4 ) 

and so 7Te,_l is subject to a two- fold degeneracy in magnitude, depending on whether the 
event impact parameters u register the lens passing on the same (— ) or opposite (+) sides 
of the lens as seen from the Earth and satellite. And 7r Ei _|_ is subject to an additional two-fold 
ambiguity in sign, depending on the relative orientations of these separations. See Gould 
(1995) for a fuller discussion of this degeneracy. 

Although no space-based parallaxes have ever been measured, this degeneracy and the 
conditions under which it can be broken have been the subject of a great deal of theoretical 
work (Gould 1995; Boutreux & Gould 1996; Gaudi & Gould 1997; Gould & Salim 1999). By 
contrast, despite the fact that more than a dozen microlens parallaxes have been measured 
using the Earth's motion (Alcock et al. 1995, 2001; Mao 1999; Soszyhski et al. 2001; Bond et 
al. 2001; Smith, Mao & Wozniak 2002; Mao et al. 2002; An et al. 2002; Bennett et al. 2002; 
Smith, Mao & Paczyhski 2003), the existence of a discrete microlens parallax degeneracy 
for the case of an accelerated observer was not discovered until the past year (Smith et al. 
2003), and the fact that there is actually a four-fold discrete degeneracy that is analogous 
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to equation (3) was only just recently recognized (Gould 2004). 

One reason for this late development is that the discrete degeneracy is typically broken 
for long events, f2 e t E > 1, where f2 e = 2n yr _1 . However, this cannot be the full explanation; 
Smith et al. (2003) showed that some previously analyzed events had unrecognized discrete 
degeneracies, and the first analysis of a four-fold degeneracy came 10 years after the event 
was discovered. Another reason, perhaps, is that while the four-fold degeneracy assumes a 
relatively simple form in the geocentric frame (Gould 2004), all previous parallax analyses 
were carried out in the heliocentric frame, in which it is much more difficult to recognize the 
physical symmetries that underlie this degeneracy. 

Whatever the exact reason for their prior neglect, it is important now to gain a better 
understanding of discrete parallax degeneracies, which are expected to be significant primar- 
ily in events that have good data and that lie below the threshold f2 e t E < 1. For cases in 
which this degeneracy can be resolved by an independent measurement of the lens-source 
relative proper motion (e.g., Alcock et al. 2001), it will be possible to measure the lens mass 
and distance. Moreover, Gould (1999) showed that the degeneracies in space-based and 
Earth-based parallax measurements are complementary, so that the full microlens parallax 
can sometimes be recovered by combining these even when each measurement is separately 
degenerate. This will have practical applications to the SIM mission, which has a major 
microlensing component. Finally, the parallax effects in the relatively short events that can 
give rise to a discrete degeneracy can also be confused with some signatures of planetary 
companions to the lens (Gaudi et al. 2002). An understanding of this degeneracy is therefore 
required for a proper analysis of microlensing planet searches. 

There are strong, although not perfect, analogies between the space-based discrete par- 
allax degeneracy first analyzed by Refsdal (1966) and the Earth-based degeneracy. In both 
cases, there is a preferred (||) axis, and the parallax component along this axis, 7Tei| is 
hardly affected by the degeneracy. For the space-based case, this axis is the Earth-satellite 
separation vector, while for the Earth-based case it is the direction of the Sun's apparent 
acceleration (in the geocentric frame). In both cases, there is a two-fold degeneracy that 
basically does not affect the magnitude of 7r E . For the space-based case, this is due to the 
first "±" in equation (4) and corresponds to a parity flip of the entire event, with the lens 
reversing the side it passes the source as seen from both the Earth and the satellite. For 
the Earth-based case, this two-fold degeneracy takes uq — > — uq and so also corresponds to 
the lens reversing the side it passes the source. Finally, in both cases, there is an additional 
two- fold degeneracy that does affect the magnitude of 7r E . For the space-based case, this is 
due to the second "±" in equation (4) and corresponds to the lens passing on the same versus 
the opposite side of the source as seen from the Earth and satellite. For the Earth-based 



- 5 - 



case, this degeneracy can be expressed in terms of the so-called "jerk parallax" , 

= 3 *V (5) 
where ex. is the (2-dimensional) acceleration of the Sun in the Earth frame projected onto 
plane of the sky and divided by an AU, and where j = dct/dt is the jerk. The degeneracy 
arises because event geometries with the same fourth-order Taylor coefficient of the lens- 
source separation, 

a 2 . o , 1 fi 2 

C 4 = — (VT E + 7Zj • 7T E ) + — — M 7Te,±, (6) 

will, to this order, generate the same lightcurve. Hence, in the limit \u \ <C 1, if 7Te is one 
solution, then 7v' E , with 

^kw = nE >\b n k± = + n jt ±), (7) 

is also a solution. 

While theoretical investigations are essential for understanding the jerk-parallax de- 
generacy, it is also important to analyze individual degenerate events. In particular, Gould 
(2004) predicted that the phenomenology of this degeneracy should be much richer for events 
observed toward the Galactic bulge than for those seen toward the Large Magellanic Cloud 
(LMC) because the bulge lies near the ecliptic while the LMC lies near the ecliptic pole. 
The jerk-parallax can be reexpressed in terms of a projected velocity, 

AU 3a 3 AU 

Vj = = r , (8) 

n jt± t E 4 a x j 

which (approximating the Earth's orbit as circular) can be evaluated, 

3 (cos 2 ip sin 2 (3 + sin 2 ip) 3/2 

4 sin/3 

Here f3 is the ecliptic latitude, ip is the phase of the Earth's orbit relative to opposition, and 
= 30kms _1 is the speed of the Earth. (Note that o: x j is a signed scalar.) Toward the 
LMC (near the ecliptic pole), Vj ~ 22kms~ 1 , independent of the time of year. However, 
toward the bulge (near the ecliptic), Vj varies from a factor sin 2 (3 below this value to | esc f3\ 
above it. Since Vj characterizes the jerk-parallax degeneracy through equations (7) and 
(8), the bulge's much broader range of values for this parameter implies a broader range of 
phenomena. 

Here we analyze MOA 2003-BLG-37, the first bulge microlensing event whose parallax 
solution has four distinct minima. We show that the analytic formalism developed by Gould 
(2004) to describe the jerk parallax applies quite well to this event. However, the morphology 
of the x 2 contours of this event is very different from that of MACHO-LMC-5, the only other 
event for which the jerk-parallax degeneracy has been detected. 



o ituo w mil u t oni w i . . 

v i = ~- znns %■ ( 9 ) 
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2. Data 

MOA 2003-BLG-37 [(RA,DEC) = (18:12:2.24,-29:01:01.3), (l,b) = (3.152,-5.656)] was 
alerted by the Microlensing Observations for Astrophysics (MOA) collaboration on 2003 
June 24 as part of its ongoing alert program (Bond et al. 2001) when the lens was about 1.5 
Einstein radii from the source. MOA obtained a total of 626 observations in I band using 
the 0.6 m Boiler & Chivens telescope at Mt. John University Observatory in New Zealand. 
Of these, 359 observations were during the 2003 season and the remaining 269 were during 
3 previous seasons. The data were reduced using on-line image subtraction. 

The Microlensing Follow Up Network (//FUN, Yoo et al. 2004) commenced observa- 
tions approximately two weeks before the peak, which occurred on HJD' =HJD— 2450000 = 
2880.68 (UT 2003 August 29.18), using the ANDICAM camera on the 1.3m SMARTS (former 
2MASS) telescope at Cerro Tololo Interamerican Observatory at La Serena, Chile. //FUN 
obtained 116 observations in / and 13 in V. The photometry was carried out using DoPHOT 
(Schechter, Mateo & Saha 1993). ANDICAM is equipped with an optical/IR dichroic beam 
splitter, allowing H band observations to be taken simultaneously with each optical image. 
However, the H band data show much larger scatter than the optical data, and so their 
incorporation does not significantly affect the solution nor reduce the errors. Hence, we 
exclude them. 

After the event peaked, we realized that it might have a measurable parallax, and we 
therefore made special efforts to observe the event regularly until late in the season so as to 
obtain good coverage of its falling wing. Since the Sun's acceleration projected on the plane 
of the sky reaches a maximum on about Sept 24 (roughly 0.6 Einstein radii after peak), the 
falling wing is the most sensitive to parallax. Observations continued several times per week 
until the end of October at both observatories. Reported errors at both observatories were 
typically 0.01 to 0.02 mag near baseline and smaller near the peak. 

Of the 626 MOA points, four were found to be mild (~ 3.5 a) outliers and were elim- 
inated from the analysis. The MOA errors were renormalized upward by a factor 1.47 in 
order to make their \ 2 P er degree of freedom (dof) consistent with unity. The //FUN / and 
V data did not contain any outliers greater than 3 a, and their x 2 /dof were both consistent 
with unity, so no adjustments were applied. 

3. Event Characteristics 

The lightcurve was initially fit to a point-source point-lens (PSPL) Paczyhski (1986) 
model, the parameters for which are shown in Table 1. It is a relatively high- magnification 
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event, with maximum magnification A max ~ 20. The data generally cover the lightcurve 
quite well. See Figure 1. Figure 2 is an instrumental color- magnitude diagram (CMD) of 
the field, which is based on yuFUN V and / data, that has been translated so that the 
centroid of the bulge clump (marked by a cross) lies at its known position (Yoo et al. 
2004), [Iq, (V — I)o]ciump = (14.32, 1.00). For stars that lie behind the disk reddening screen, 
the diagram therefore shows their dereddened color and magnitude, while foreground stars 
that lie within the dust distribution (e.g., the "reddening sequence" toward the upper left) 
appear in the diagram as brighter and bluer than they actually are. The dereddened color 
and magnitude of the microlensed source (when unmagnified) are shown as an open circle. 
These values are determined from the fits for the deblended source fluxes f s in V and / 
from /iFUN, i.e., the same data set used to construct the CMD. Using jackknife, we find 
that the statistical errors from centroiding the clump are 0.02 mag in each axis. The entire 
procedure assumes a uniform reddening over the 6' x 6' field. In fact, we find that when 
we consider stars within 100" of the source star (and so shrink the sample size by 75%), 
the clump centroid shifts by only 0.02 mag in color and 0.05 mag in brightness, which are 
both consistent with statistical fluctuations. That is, there is no evidence for differential 
reddening on the relevant scales. 

Based on its position in the CMD, the source could be a clump giant in the bulge but 
it could equally well be a disk first-ascent giant that lies ~ 1 mag in the foreground. We 
determine the angular size of the source from its position on the instrumental CMD using 
the standard approach as summarized by Yoo et al. (2004). Very briefly, we use standard 
color-color relations to convert the observed (V — I) = 0.91 to (V — K) = 2.06. Then 
using the dereddened V = 15.27 and the van Belle (1999) empirical color /surf ace-brightness 
relation, we find, 

0* = 5.8 ± 0.5 //as. (10) 

The error is almost entirely due to the 8.7% scatter in the van Belle (1999) relation. The 
next largest source of error comes from the 0.07 mag uncertainty in the fit to the source 
flux, which contributes about 3% in quadrature to 6*. Other sources, such as the 0.004 
mag error in the fit color and the 0.02 mag errors in the color and magnitude of the clump 
centroid, contribute still less. The total error is 9.5%. (Note that this entire procedure rests 
on instrumental magnitudes, so absolute calibration of the photometry is not necessary.) 

We consider finite-source effects, holding p = 9*/9e at various values while fitting for 
the other parameters. We find a minimum x 2 — 751.37 at p = 0.023, but since this is only 
Ax 2 = 1-7 lower than at p = (see Table 1), we do not consider this to be a significant 
detection of finite-source effects. We are, however, able to put an upper limit on p, and so 
a lower limit on 6>e = 0*/ p and hence on the lens-source relative proper motion p = 0e/^e- 
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We find, 

p < 0.034, 6> E > 170,uas, p, > 1.5masyr _1 , (11) 

at the 2 a level. (At 3 a, p < 0.44.) These constraints are of marginal interest. The typical 
proper motion expected for bulge-bulge lensing events is about 5 masyr^ 1 , and for disk-bulge 
events it is about 8masyr _1 . Hence, only the low-/i part of parameter space is eliminated. 
Similarly, from equation (1) for 9 E , we obtain, 

M 7T rP i . , 

> 0.15. (12) 



M 25/xas 

Hence, only the very low-mass lenses with very small lens-source separations are eliminated. 
Microlensing flux observations F(t) are always fit to the form, 

F(t) = f s A(t) + f b , (13) 

where A(t) is the model magnification as a function of time, f s is the baseline flux of the 
lensed source, and ft, is the background light within the point spread function (PSF) that 
does not participate in the microlensing event. A somewhat unsettling feature of the PSPL 
model is that it has negative blending f b in both I and V, which is detected at the 6 a level 
and so cannot be due to ordinary statistical fluctuations. Logically, there are only three 
other possible explanations: a failure of the model to take account of some essential physics 
affecting the lightcurve, systematic errors in the photometry, or a "negative source of light" 
within the PSF. 

The last explanation is not as absurd as it might appear. The dense star fields of the 
Galactic bulge contain a mottled background of upper-main-sequence and turnoff stars that 
lie below the detection limit. If the source happens to lie in a "hole" in this background, the 
fit to the blended flux will be negative. By comparing the negative blend fluxes, fb-, to the 
source fluxes f s , and referencing this comparison to the CMD, we find that the "hole" has a 
dereddened color and magnitude of (J, V — I)o,"hoic" ~ (17.1,0.5), which is about the right 
color but somewhat too bright to be a turnoff star. The position of this "hole" is shown as 
a filled triangle in Figure 2. 

While systematic errors are possible, it seems unlikely that they would afflict the in- 
dependently measured V and / photometry by similar amounts. The remaining possible 
explanation is an inadequate model. Since blending is even in time (about the event peak), 
it can easily be confused with other even effects. The two most obvious possibilities are 
finite-source effects and parallax. The above mentioned 1.3 a "detection" of finite-source 
effects does indeed tend to reduce the negative blending, but by a negligible amount. The 
7Te,|| component of 7r E gives rise to a third-order term, which is asymmetric in time and hence 
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cannot be confused with blending. However, the tte,± component gives rise to a fourth-order 
term, which can be confused with blending (Smith et al. 2003). 

4. Parallax Fits 

The lightcurve shows a very minor asymmetry, being slightly brighter than the model 
on the rising side and slightly fainter on the falling side. See Figure 1. This, together with 
the negative-blending problem described in § 3, lead us to fit the lightcurve for parallax. We 
employ the geocentric formalism of Gould (2004). We minimize x 2 using Newton's method, 
starting with the no-parallax solution as a seed. We find a solution whose x 2 is lower by 
approximately 32, which we regard as significant. See Table 1. To find the corresponding, 
uq < solution predicted by Smith et al. (2003), we adopt a seed solution with uq — > — uq 
and with all other parameters the same as the uq > solution. We find that x 2 is minimized 
extremely close to this seed. The difference in x 2 of these two solutions is less than 1. See 
Table 1. 



4.1. Jerk-Parallax Predictions 

To navigate the jerk-parallax degeneracy, we employ the formalism of Gould (2004). 
We first evaluate the acceleration and jerk of the Sun projected onto the plane of the sky at 
time t - We use the method of An et al. (2002) to find the path of the Sun projected on the 
plane of the sky and take numerical derivatives to obtain in (north,east) coordinates, 

a = (0.060, 0.882)fi|, j = (-0.075, 0.441)fi|. (14) 

Hence, the solar acceleration points almost due east, with a position angle of 86° (north 
through east). Since en determines the parallel direction, the southern u > solution for 
7Te listed in Table 1 can be represented in the (||, _L) coordinate system as, 

tte = (7Te.ii, 7t e ,±) = (-0.049, 0.568), (15) 

This explains why the parallax signature in Figure 1 is so weak. The component of 7Te 
parallel to the solar acceleration, which is what gives rise to the parallax asymmetry, is 
extremely small. This asymmetry arises from the third-order term in the Taylor expansion 
of the lens-source separation (Gould 2004). Since this is the lowest order affected by parallax, 
it is also the easiest to notice in the lightcurve (Gould, Miralda-Escude & Bahcall 1994). 
The perpendicular component first enters at fourth order (Gould 2004) and therefore is much 
harder to recognize. 
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From equations (5), (8), and (14), we find 

n j>± = 0.274, or iJj = 147kms _1 . (16) 

Combining equations (7), (15), and (16), and then switching to (north,east) coordinates, we 
predict that the degenerate solution will lie at, 

«north, <cast) - (0-843, -0.009) (predicted). (17) 

We begin with this as a seed, and find that \ 2 is minimized at (ir' E north , 7r Ejeast ) = (0.788, —0.023), 
which is quite close to the predicted value. See Table I. 

We perform a similar exercise for the u < solutions. We predict the degener- 
ate solution will lie at (^.north) ^E.east) — ► (0.863,-0.017) whereas it is actually found at 
(^E.north; ^E.east) = (0-783, —0.012). We conclude that the formalism predicts the location of 
the jerk-parallax degenerate solutions quite well, at least in this case. 

Some geometrical insight into this procedure can be gained from Figure 3, below. Ac- 
cording to equation (7), the midpoint of the two solutions should lie at 

^,^ = ^^ = (^,11,-^/2). (18) 

The first component says that both solutions (as well as their midpoint) should lie along a 
line defined by the Sun's apparent acceleration (labeled 7Te,|| in the figure). This is indeed 
the case. The second component says that the midpoint should lie along this axis and 
— 7r j,±/2 = —0.137 from the "perpendicular axis" (which is parallel to the arrow labeled 7r Ei _i_ 
and passes through the origin). This is also very nearly the case. 



4.2. Timescale Predictions 

The jerk-parallax formalism predicts that of the remaining parameters (t , Uo, ^e, fs, fb), 
all but u and t E should be almost the same for all four solutions. Moreover, uq should only 
change sign, while t E should differ among solutions by, 

At E ~ -^(vr Ei± AM + M A7T Ei± )t E . (19) 

Table 1 shows that these predictions are also borne out. The parameters to, \v>o\, and the 
three sets of (f s , fb) are all virtually unchanged among the four solutions. Equation (19) 
predicts At E = (0.60, 0.89, 0.46) days when comparing, respectively, the two southern, two 
northern, and corresponding southern and northern solutions. The actual differences from 
Table 1 are (0.62,0.83,0.70) days for the southern and northern pairs, and for the average of 
the two south/north comparisons, respectively. 
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4.3. 



Parallax Contours 



To further explore the degeneracy, we calculate x 2 f° r a grid of parallax values and show 
the resulting likelihood contours in Figure 3. There are three features to be noted. First, 
there are two local minima, which as discussed above are related to each other almost exactly 
as predicted by equation (7). Second, the northern minimum is higher by Ax 2 = 7.5 and so 
is formally excluded at the 2.6 a level. See also Table I. Third, each minimum is elongated 
along the 7Te,± axis and, indeed, is embedded in a larger continuous degeneracy also along 
this axis. That is, at the 4 a level, there is a continuous set of acceptable solutions that 
extends almost from vr E north = —1 to vr E north = 1, but whose width is less than 0.005. The 
two discrete solutions form islands within this linear degeneracy. 



These three features should be contrasted to the case of MACHO-LMC-5, for which 
the analogous contour diagram is shown in figure 3 of Gould (2004). For that event also, 
the jerk-parallax formalism accurately predicted (or rather postdicted) the location of the 
second minimum, and consequently the separation vector between the two solutions was 
perpendicular to the Sun-acceleration vector. However, in other respects the two cases are 
quite different. First, the two MACHO-LMC-5 minima had nearly identical x 2 's- Second, 
these two minima were not severely elongated. Third, they did not lie within the continuous 
degeneracy. 

In fact, Gould (2004) argued that there should be a hierarchy of microlens parallax 
degeneracies. When the data were adequate to solve for the source-lens separation only to 
second order in time, there would be no parallax information. When they were adequate to 
solve the separation to third order, there would be a continuous (linear) degeneracy. When 
they were adequate to solve to fourth order, the continuous degeneracy would be broken, but 
a discrete degeneracy would remain. Finally, if higher order terms could be measured, the 
degeneracy would be completely lifted. This picture would seem to argue that the discrete 
degeneracy should be embedded in the continuous degeneracy as it is in the present case of 
MOA 2003-BLG-37. It is curious therefore that for MACHO-LMC-5, the degeneracy is not 
embedded, even though it is more severe (in the sense the two minima have almost identical 



5. 



Discussion 



5.1. Comparison to MACHO-LMC-5 
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5.2. Negative Blending 

After including parallax in the fit, negative blending remains in both V and /, but it is 
less severe. The dereddened color and magnitude of the "hole" is, 

Vhoic" = 17.7 ± 0.2, (V - /) ,"hoie» = 0.2 ± 0.3, (20) 

where the errors are approximately independent. The hole still has the right color (within 
errors) to be due to fluctuations in the turnoff-star background, and its amplitude is substan- 
tially smaller and so is consistent with an upper main-sequence or turnoff star. See Figure 2 
(open triangle). 

5.3. Xallarap 

As Smith et al. (2003) point out, any lightcurve that can be fit by parallax can also be 
fit by xallarap (accelerated motion of the source). For long events with excellent data, the 
parallax parameters are so well determined that the putative xallarap solution can succeed 
in mimicking parallax only by having the source mirror the projected motion of the Earth 
to extremely high precision. In such cases, one may apply Occam's razor to argue that the 
xallarap solution is unlikely. For some other events, there may be external information that 
tends to corroborate the parallax interpretation, even if the parallax solution itself is not 
tightly constrained. For example, the lens of MACHO-LMC-5 was directly imaged which 
directly confirmed the direction of motion (Alcock et al. 2001) as well as the lens mass and 
distance (Gould 2004) derived from the parallax solution. 

In the present case, the only auxiliary information comes from the constraint #e > 
170yuas, which was derived from the lack of finite-source effects. See equation (11). If the 
source is assumed to be at D s = 8kpc, this implies that the Einstein radius projected on 
the source plane is f E = D s 0e > 1.36 AU. To mimic the effect of the Earth's accelera- 
tion, the source must have a projected acceleration a S)Pro j such that a SiPr oj/^E = ae.proj/^E, 
where a eiPro j = aAU is the projected acceleration of the Earth. That is, a SjPro j = K E ar E = 
TTEaD s 9 E > 0.6a e , where a® = AUfi| is the (full) acceleration of the Earth. Such a source 
acceleration is of course possible. However, if such an acceleration were induced by a compan- 
ion of mass m, then the velocity semi-amplitude would be at least v > 23 (m/M G ) 1//2 kms _1 . 
Hence, unless the companion were substellar or the system were seen pole on, it would be 
possible to detect such motion by radial velocity measurements. 
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6. Conclusions 

We have shown that MOA 2003-BLG-37 suffers from a four-fold jerk-parallax degener- 
acy. The two pairs of solutions have micro lens parallaxes of 7Te = 0.57 and 7r E = 0.78. These 
correspond to projected Einstein radii of f E = AU/7r E = 1.76 AU and f' E = 1.28 AU. The 
first solution is favored at the 2.6a level. We were able to place only lower limits on the 
angular Einstein radius, 6>e > 170 //as. Combining the favored parallax measurement with 
this limit yield lower bounds on the mass and lens-source relative parallax of 

M > 0.037 M Q , 7r rcl > 100 //as, (21) 

neither of which probes a really interesting regime. 

MOA 2003-BLG-37 is only the second event for which four separate minima have been 
found in x 2 as a function of the vector microlens parallax 7r E . The other was MACHO- 
LMC-5. For both events, the jerk-parallax formalism of Gould (2004) accurately predicts 
the position of the second pair of solutions given that the first pair has been found. This 
gives hope that the additional solution will easily be found in the future events that contain 
them. 

On the other hand, there are substantial differences between the morphologies of the 
X 2 contour diagrams for these two events in terms of the shapes and relative depths of the 
minima as well as the degree to which the minima are embedded in the lower-order contin- 
uous degeneracy. Further understanding of these differences is likely to require additional 
theoretical work as well as the study of additional events that display the jerk-parallax de- 
generacy. In particular, it would be useful to check all events with measured parallaxes to 
determine which of these are affected by this degeneracy. The black-hole candidates of Mao 
et al. (2002) and Bennett et al. (2002) would be especially interesting in this regard. 
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Fig. 1.— Lightcurve of MOA 2003-BLG-37 showing data from MOA / (black), /iFUN I 
(red), and //FUN V (green). The fit does not allow for parallax. The residuals show a very 
slight asymmetry, with the data brighter than the model on the rising wing and fainter on 
the falling wing. The MOA data are shown as recorded, and each /xFUN dataset has been 
linearly rescaled to be consistent with the MOA system. The MOA scale is not rigorously 
calibrated but is accurate to a few tenths of a magnitude. 
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Fig. 2. — Instrumental color-magnitude diagram of the 6' x 6' field centered on MOA 2003- 
BLG-37. The diagram has been translated so that the centroid of the bulge clump giants 
(cross) is at its known position [J , (V — 7)o] c iump = (1-00, 14.32). The open circle shows the 
color and magnitude of the (deblended) unmagnified source as determined from the fit to the 
lightcurve. Using this diagram and the standard technique described in Yoo et al. (2004), 
we find that the source has an angular radius of 6* = 5.8 ± 0.5 //as. The solid triangle shows 
the color and magnitude of the "hole" in the mottled background of bulge turnoff stars that 
is required to account for the negative blending in the no-parallax model. When parallax is 
included, the magnitude of this "hole" (open triangle) shrinks, and is more in line with what 
is expected. 
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Fig. 3. — Likelihood contours for the u < solutions in the microlens parallax (7r E ) plane 
shown at A% 2 = 1, 4, 9, 16, 25, 36, and 49 relative to the minimum. The magnitude 
of 7Te gives the size of the Earth's orbit relative to the Einstein radius, and its direction 
is that of the lens-source relative motion. There are two local minima, with the northern 
minimum higher by A% 2 = 7.5 and hence surrounded by the 3 a contour. The two minima 
are embedded in contours of the continuous line-like degeneracy, which runs perpendicular 
to the direction of the Earth-Sun acceleration vector (arrows) and which reflects the parallax 
asymmetry. Note that since 7Te is dimensionless, the axes have no unit labels. This diagram 
should be contrasted with figure 3 of Gould (2004) as described in § 5.1 of this paper. 



Table 1. Jerk-Parallax Degenerate Solutions of MOA-2003-BLG-37 



to 






tte,jv 




F s ,MOA-I 


Fb,MOA-I 


-Ps.uFUN-/ 


-Ffc.uFUN-/ 


F s ,nFVN-V 


Fb,uFVN-V 


x 2 


days 




days 






(1) 


(1) 


(2) 


(2) 


(3) 


(3) 




2880.6728 


0.051488 


42.5576 


0.0000 


0.0000 


9.2533 


-0.7179 


29.1579 


-2.2273 


12.5789 


-1.4197 


753.093 


2880.6791 


-0.050506 


43.5304 


-0.5882 


0.0818 


9.0804 


-0.5424 


28.5826 


-1.2718 


12.3363 


-1.0411 


721.054 


2880.6797 


0.050509 


42.9087 


-0.5629 


0.0872 


9.0809 


-0.5423 


28.5853 


-1.2748 


12.3370 


-1.0394 


720.719 


2880.6787 


-0.050756 


42.6415 


0.7829 


-0.0124 


9.1258 


-0.5932 


28.7126 


-1.4132 


12.3949 


-1.1167 


728.594 


2880.6778 


0.050849 


43.4089 


0.7881 


-0.0231 


9.1426 


-0.6111 


28.7639 


-1.4702 


12.4168 


-1.1392 


730.522 


M 



I 



Note. — 1 = In units of / = 18, where / is the (approximately calibrated) MOA apparent mag shown in Figure 1. 
Note. — 2 = In units of Iq = 18, where Iq is the dereddened apparent mag used in Figure 2. 
Note. — 3 = In units of Vq = 18, where Vq is the dereddened apparent mag used in Figure 2. 



